EPJ manuscript No. 

(will be inserted by the editor) 



The pion-nucleon scattering lengths from pionic deuterium 



in 
o 
o 

> 
O 

o 

cn 

> 
cn 

o 

^ 

o 
o 



Ulf-G. Meifiner^'^ Udit Raha\ and Akaki Rusetsky^'^ 

^ Universitat Bonn, Helmholtz-Institut fur Strahlen- und Kernphysik (Theorie), 

Nufiallee 14-16, D-53115 Bonn, Germany 
^ Forschungszentrum Jiilich, Institut fiir Kenphysik (Theorie), D-52425 Jiilich, Germany 

On leave of absence from: High Energy Physics Institute, Tbilisi State University, 

University St. 9, 380086 Tbilisi, Georgia 

Received: date / Revised version: date 

Abstract. We use the framework of effective field theories to discuss the determination of the S-wave 
■kN scattering lengths from the recent high-precision measurements of pionic deuterium observables. The 
theoretical analysis proceeds in several steps. Initially, the precise value of the pion-deuteron scattering 
length a-nd is extracted from the data. Next, a-^d is related to the S-wave nN scattering lengths a+ and 
a_. We discuss the use of this information for constraining the values of these scattering lengths in the full 
analysis, which also includes the input from the pionic hydrogen energy shift and width measurements, 
and throughly investigate the accuracy limits for this procedure. In this paper, we also give a detailed 
comparison to other effective field theory approaches, as well as with the earlier work on the subject, 
carried out within the potential model and multiple scattering framework. 
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1 Introduction 

The Pionic Hydrogen collaboration at PSI [1,2] has per- 
formed high precision measurements of the strong interac- 
tion shift Els and width Fis of the Is state of pionic deu- 
terium from the 3p— Is X-ray transition. The (complex) 
pion-deuteron scattering length was extracted from these 
measurements with the use of the leading-order Deser for- 
mula [3] 



. - „ -a^d, (1.1) 

2 rs 

where Eis — \ a^f^d is the Coulomb binding energy in 
the Is state, and rs = {a^d)~^ denotes the Bohr ra- 
dius (in these formulae, Hd stands for the reduced mass 
of the wd system). The most recent measurement of the 
pion-deuteron scattering length by the Pionic Hydrogen 
collaboration at PSI [2] yields 



-0.0261 (±0.0005) + i 0.0063 (±0.0007) M, 



(1.2) 

Performing experiments to determine aT^d is usually jus- 
tified by the possibility of extracting independent infor- 
mation about the ttN 5- wave isoscalar (a+) and isovector 
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(a_) scattering lengths. What makes this enterprise par- 
ticularly interesting is the fact that in the multiple scat- 
tering theory a^^ = const • a+ + correction terms. If one 
could accurately evaluate the higher-order terms in this 
expression, then a precise measurement of a^^ would en- 
able one to constrain the value of a+, which is in general 
a rather delicate task. The reason for this is that, since 
a+ is much smaller than the isospin-odd scattering length 
a_, a very high accuracy is needed in order to determine 
a+ from the measurements of the linear combinations of 
a+ and a_. This is exactly the case in the experiments 
measuring the pionic hydrogen energy shift and width, 
which enable one to determine the combinations (a++a_) 
and a_ , respectively. Thus, the measurement of a-^d con- 
tributes a complementary piece of information about the 
scattering lengths, which can be used in the complex theo- 
retical analysis finally aimed at the determination of both 
a+ and a_ [4]. We wish to mention here that these scat- 
tering lengths are quantities of fundamental importance 
in low-energy hadronic physics, since they test the QCD 
symmetries and the exact pattern of the chiral symmetry 
breaking. Moreover, since the knowledge of these scatter- 
ing lengths places a constraint on the ttN interactions at 
low energy, it also affects our understanding of more com- 
plicated systems where t:N interaction serves as an input, 
e.g. NN interaction, 7r-nucleus scattering, three-nucleon 
forces, etc. 

Expressing the 7rd scattering length in terms of the pa- 
rameters characterizing the underlying pion-nucleon and 
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nuclcon-nuclcon dynamics is one of the classical problems 
in conventional nuclear physics based on the potential 
scattering formalism, see, e.g. [5-13]. Note, however, that 
the experiments on pionic deuterium will be used to ex- 
tract the ttN scattering lengths in QCD and not in any 
potential model. In other words, using the latter in order 
to establish the relation between nd and ttN scattering 
lengths introduces a theoretical error in the analysis of 
the experimental data, whose magnitude is very hard to 
control. 

In recent years, the problem of a very low energy pion- 
deutcron scattering has been studied within the frame- 
work of effective field theories (EFTs). The method orig- 
inates from the seminal paper of Weinberg [14], where 
chiral Lagrangians have been systematically applied for 
the description of interactions of pions with nuclei. In this 
paper, by using the chiral Lagrangian, one calculates the 
set of diagrams contributing to the "irreducible transition 
kernel" for the pion scattering on two nucleons, and the re- 
sult is then sandwiched between "realistic" dcuteron wave 
functions in order to evaluate the scattering amplitude of 
the pion on the deuteron. In actual calculations carried out 
in Ref. [14], the phenomenological deuteron wave function 
for the Bonn potential model has been used. We wish to 
note that this last step, in general, can be justified within 
the framework of the effective field theories, only if the 
particular process which one is going to describe is domi- 
nated by the long-range mechanisms, e.g. by one-pion ex- 
change. On the other hand, when the calculations within 
such a "hybrid" approach arc pursued in higher orders 
in Chiral Perturbation theory (ChPT), the kernels grow 
faster with a large momenta and probe shorter distances. 
Moreover, this short-distance behavior is not necessarily 
correlated to that of the phenomenological wave function. 
From this we conclude that in order to obtain a systematic 
description of the pion-deuteron system, it is preferable to 
use the deuteron wave functions and the transition ker- 
nels, evaluated within the same field-theoretical setting - 
in this case, no specific conjecture about the dominance 
in the unobservablc quantities, like the kernel or the wave 
function is needed. For the latest work within the hybrid 
approach, sec e.g. [15, 16]. 

Further development of the approach based on chi- 
ral Lagrangians (sec, e.g. [15-20]), has followed different 
paths. In the paper [18] one has used the framework with 
perturbative pions, whereas the authors of Ref. [19] make 
use of the so-called Heavy Pion EFT (HP EFT) with 
the dibaryon field. The latter approach is quite close to 
the one used in the present work. The technique used 
in Refs. [18, 19] has the advantage that one may easily 
construct the deuteron wave function in a closed form, 
since the lowest-order nucleon-nucleon interactions are de- 
scribed by a contact four-nucleon vertex. The central prob- 
lem in both the papers is related to the calculation of one 
particular diagram, describing double scattering of pions. 
These calculations lead to a very strong scale dependence 
near /x ~ M^r - a natural choice of the scale parameter 
in this sort of the effective theories. Since, on the other 
hand, this dependence must be canceled by a contribu- 



tion from the low-energy constant (LEG), which describes 
pointlike interactions of four nucleons and two pions, we 
easily conclude that the magnitude of this LEG can not be 
small. In the absence of any information about the actual 
value of this LEG apart from naive dimensional order-of- 
magnitude estimates based on the naturalness arguments, 
one may finally conclude that the theoretical uncertainty 
in the relation of nd and ttN scattering lengths should be 
very large. 

On the other hand, the results obtained in Ref. [20] 
seem not to be in agreement with those of Refs. [18, 19]. 
The method which is used in Ref. [20] , is a systematic ex- 
tension and elaboration of Weinberg's original proposal, 
where both the transition kernel and the dcuteron wave 
function arc constructed in ChPT (note also, that the sys- 
tematic derivation of the unitary and the energy-indepen- 
dent potentials within this framework has been discussed 
recently in Ref. [21]). The approach uses cutoff regulariza- 
tion to deal with the potentials that are growing for a large 
three-momenta. The typical scales for the cutoff mass A 
are somewhat smaller than the hadronic scale in QCD 
^ 1 GeV (depending on the order in ChPT in which the 
calculations are carried out). The results of the calcula- 
tions are /l-dcpendcnt, which is a reminiscent of the scale- 
dependence in the dimensional regularization scheme. The 
bulk of this yl-dependence should be canceled by an analo- 
gous dependence in the LEGs, and the remainder, which is 
an artifact of the non-perturbative formalism used, should 
be of a higher order in ChPT. In Ref. [20], the cutoff de- 
pendence of the nd scattering length has been studied, 
with the LEGs assumed to vanish. In a remarkable con- 
trast with Refs. [18,19], the yl-dependcnce of the results 
in Ref. [20] turns out to be very mild, thereby conclud- 
ing that the LEGs must have a weak cutoff dependence. 
If one could interpret the cutoff dependence as an esti- 
mate of the uncertainty of the method, then the results of 
Ref. [20] would amount to a rather accurate prediction of 
the Tvd scattering length within the framework of ChPT. 

The present situation which was described above is 
unacceptable from the point of view of both the theory 
and the phenomenological analysis of the data. From the 
theoretical point of view, the calculations carried out in 
Ref. [20], clearly indicate that the diagrams in which the 
virtual pions are emitted or annihilated, arc strongly sup- 
pressed. This phenomenon originates from the infrared 
enhancement of a certain class of the diagrams in the 
Weinberg scheme, as well as the threshold suppression of 
the diagrams containing pseudoscalar vertices. In order to 
accommodate the above feature in the theory with non- 
perturbative pions, in Ref. [20] a novel counting, inspired 
by the HP EFT, has been enforced on top of the conven- 
tional ChPT Lagrangian. Stated differently, this means 
that simpler effective theories, which were used in Refs. 
[18,19], are physically adequate for the problem consid- 
ered. How can it then be that using a simpler theory, we 
get an answer which contradicts the answer obtained in 
Weinberg's framework [20], the very approach one starts 
from? From the point of view of phenomenology, the ex- 
isting confficting predictions, on the one hand, do not en- 
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courage the experimentalist's efforts to measure the pion- 
deuteron scattering length to a better accuracy and on the 
other hand, suggest that the values of the ttN scattering 
lengths extracted from the analysis of the nd data should 
be taken with a grain of salt. 

The aim of the present paper is to perform a thor- 
ough investigation of pion-deuteron scattering at thresh- 
old within the framework of low-energy effective theories. 
In particular, we plan to clearly establish the limits of 
accuracy for extracting 'kN scattering lengths from the 
measured ttc? scattering lengths. We also perform a de- 
tailed study of the above-mentioned discrepancy between 
the results obtained within the HP EFT and in the Wein- 
berg approach. Moreover, the investigation of this subtle 
question, in our opinion, is by itself very informative and 
sheds light on many peculiar aspects of the effective field 
theories in general. 

The complex problem, which we are going to consider 
in this paper, naturally falls into several sub-problems, 
which are characterized by distinct momentum scales. 
Consequently, instead of trying to describe everything at 
once, it is convenient to construct a tower of effective field 
theories, matched one to another, each designed for one 
particular momentum scale. 

i) At the momentum scales afSd — 1 MeV, the charged 
pion and the deuteron form an atom, whose observ- 
ables are measured by the experiment. The character- 
istic distances in such an atom - hundreds of fm - 
are much larger than the deuteron size, and the bind- 
ing energy in the ground state, which almost coincides 
with the Coulomb binding energy Eis — ^ cP' — 
3.5 • 10'^ eV, is much smaller, than the binding en- 
ergy of the deuteron e = 2.22457 MeV. Stated differ- 
ently, at these energies the deuteron can not be re- 
solved as a composite particle, and the effective the- 
ory, which describes the atom, contains the deuteron 
field (not the nuclcon fields) as an elementary degree 
of freedom. The hard momentum scale in this effec- 
tive theory is given by the average value of the three- 
momentum of the nucleons bound within the deuteron, 
7 = ^Jern ~ 45 MeV, where m stands for the nucleon 
mass. The expansion parameter in this theory is given 
by the ratio of the scales ay^dll = 0(a). The out- 
put from the calculations within this effective theory 
is relation which connects the measured energy shift of 
the 7rd bound state to the 7rd scattering amplitude at 
threshold in the next-to-leading order in isospin break- 
ing. In its turn, the latter at the leading order in isospin 
breaking coincides with the 7rd scattering length a-i^d- 

ii) Extracting the scattering length from the pionic deu- 
teron one next has to find the relation of this quantity 
to the ttA^ 5- wave scattering lengths a+ and a_. In 
order to achieve this goal, we have to construct an- 
other effective field theory, in which the independent 
degrees of freedom are the pion and the nucleon fields, 
whereas the deuteron emerges as a bound state of the 
proton and the neutron. The characteristic momen- 
tum scale in this theory is defined by the binding mo- 
mentum 7. Furthermore, a careful analysis of the re- 



sults of Ref. [20] provides us with an important clue: 
The processes, in which the virtual creation and anni- 
hilation of pions takes place, are suppressed as com- 
pared to the processes where this does not occur (al- 
though both processes may formally have the same chi- 
ral order). Note that these processes naturally come to- 
gether in the conventional relativistic QFT. This fact 
clearly indicates that the most economic way to de- 
scribe Trd scattering at threshold is to design an effec- 
tive field theory, in which the pion creation and annihi- 
lation processes are explicitly excluded - all vertices in 
the Lagrangian contain equal number of ingoing and 
outgoing pions and nucleons. The whole information 
about these processes is, however, not lost: it is in- 
cluded in the pertinent LECs of such an effective the- 
ory. Moreover, it is also clear that for such small en- 
ergies, one can treat kinematical relativistic factors as 
perturbations both for pions and for nucleons. 
The calculations of the deuteron properties in the 
above-described theory, which will be referred to as 
the heavy-pion effective theory hereafter, dramatically 
simplify and can be performed analytically. The out- 
put of the calculations is the quantity a.,^di expressed in 
terms of the threshold parameters of the -kN and 
scattering. The hard scale in such a theory is given 
by the pion mass M^r , and the expansion parameter is 
given by the ratio of scales ^/M^^ ~ 1/3. The matching 
to the previous effective theory is performed for the -nd 
scattering amplitude at threshold: this quantity must 
be the same in both theories. 

Note also that from now on we neglect all isospin- 
breaking effects (one could not do this at the earlier 
step, because the pionic deuterium is created predom- 
inately by Coulomb interactions.). In this approxima- 
tion, the threshold scattering amplitude coincides with 
the TTc? scattering length a-j^d- If needed, the isospin- 
brcaking effects can be turned on later, 
iii) The simplicity of the calculations in the HP EFT co- 
mes at the cost of the large size of the LECs. Since 
the hard scale of the theory is determined by M^r, 
this is also the scale that enters in the estimate of 
the size of the (unknown) LECs in the assumption 
that these LECs have the natural size (note that some 
LECs might be parametrically enhanced as compared 
to the value which is expected on the purely dimen- 
sional grounds, see below). On the other hand, if the 
calculations are done in ChPT, the natural-size LECs 
are suppressed by a higher scale 47ri^7r — 1 GeV rather 
than A/tt. Thus, the rationale for performing calcula- 
tions in the Weinberg framework can be formulated as 
follows. In these calculations, one "resolves" the dy- 
namics of the system at the scales from M^^ up to the 
scale ^ 1 GeV, which is the energy range where the in- 
teractions in the system of few pions and nucleons are 
predominately determined by (multi)-pion exchanges. 
One may then assume that the bulk contribution to 
the HP LECs comes from the momentum region be- 
tween Mjr and 1 GeV and can be expressed in terms 
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of pion loops, which arc calculated in the Weinberg 
scheme. 

If we suppose that such a scheme is realized, we ar- 
rive at the effective theory, where the characteris- 
tic momenta are of order and the hard scale 
is given by 47rF^. The expansion parameter, in the 
absence of other scale, is given by the ratio of two 
scales M^/{A'kF^). The matching to the HP EFT is 
performed for the iS-matrix element of the process 
ttNN ttNN, that determines a particular LEG 
of the HP Lagrangian. One of the objectives of the 
present paper is to find out whether doing the calcula- 
tions in the Weinberg framework and performing the 
matching to the HP EFT enables one to indeed reduce 
the uncertainty related to the choice of LECs. 

The organization of the paper follows the above-des- 
cribed scheme of "nested" effective field theories. Namely, 
in section 2 we consider the precise extraction of the nd 
scattering length from the experimental data on the pionic 
deuterium. Then, in section 3, we construct the systematic 
heavy-pion effective theory (HP EFT) in order to calculate 
the nd scattering lengths in terms of the threshold param- 
eters of ttN and NN interactions. In order to establish the 
connection to ChPT in the Weinberg scheme, in section 4 
we perform the matching of the threshold amplitudes in 
both theories. We also provide a numerical analysis and 
discuss the question of accuracy. A detailed comparison to 
the existing approaches is carried out in section 5. Finally, 
section 6 contains our conclusions. 



2 Pionic deuterium 

In the experiment at PSI [1,2], one measures the energy 
of 3p — Is X-ray transition, deducing the strong shift of 
the pionic deuterium in the Is state and the nd scatter- 
ing length from this measurement. At the first step, in 
order to obtain the strong shift, one has to subtract the 
so-called "electromagnetic shift" from the full measured 
value, where the former is calculated in the accuracy that 
matches the experimental precision. At the next step, the 
Trd scattering length should be extracted from the strong 
shift by means of the Deser-type formula (1.1). If required 
for accuracy considerations, the latter relation can also 
be generalized to include next-to-leading order isospin- 
breaking corrections. 

To the best of our knowledge, complete calculations of 
the electromagnetic shift in the pionic deuterium are not 
available in the literature, except the results contained in 
table 1 of Ref. [1] , where different contributions are given 
without a derivation. The investigations in Ref. [22] are 
not complete - as the authors themselves note, they do not 
include all isospin-breaking corrections at next-to-leading 
order. In order to have a complete and transparent field- 
theoretical treatment of the pionic deuterium problem at 
all levels, we find it appropriate here to re-derive the ex- 
pression for the full energy shift at order a*, a^{md — mu), 
and to check (at least, numerically) the results given in ta- 
ble 1 of Ref. [1]. 



The method, which will be used in our calculations, is 
analogous to the one applied recently to describe tt'^tt^ 
[23-26], TT-p [27,28], ttK [26] and R-p [29] atoms. In this 
section, we display only the final results of the calculations 
- the necessary details are provided in appendix A. The 
full binding energy in a given stationary state of the pio- 
nic deuterium depends on the principal quantum number 
n, on the orbital quantum number I and on the total an- 
gular momentum j. For a given / (except / = 0) the total 
angular momentum j takes the values j = I ~ 1,1,1 + 1. 
This splitting, which is explicitly evaluated in appendix 
A, is tiny. The following averaged value is relevant for the 
analysis of the experimental data 



Enl — 



1 



6(1 + 1/2) 



(2.1) 



Up to the next-to-leading order in isospin breaking, the 
full energy shift of the nl state can be separated in what 
is called "electromagnetic" and "strong" parts. In order to 
simplify the comparison to the existing results, the former 
is additionally split by hand in different pieces. Finally, at 
this order one obtains 



Enl 



Enl - 
Enl 



(2.2) 



In the above formula, we have chosen the same naming 
scheme as in Ref. [1]. Note that in this paper individual 
contributions are not specified explicitly, so the identifica- 
tion, which is given below and in table 1, is performed by 
analogy with the pionic hydrogen case [30]. Our explicit 
expressions are given below 



KG 



' 2n2 



1 



3aV^ 



2n 
21 + 1 



An 



8n^{Md + Mn) \l + l/2 
4n6io - 4 



3 

' 4 

- 3 



6?! 



AMdM^n^ V 1 + 1/2 

AEll^6l,^^a'^,l{{rl) + {rl)), 



(2.3) 



where Md denotes the mass of the deuteron. An explicit 
expression for the vacuum polarization contribution is gi- 
ven in Ref. [31], see Eq. (3) of that paper. 

In order to be able to compare with the existing results, 
our numerical calculations have been performed for the 
same values of the input parameters as in Ref. [1]. We take 
the deuteron binding energy to be e = 2.22457 MeV, and 
the charge radii of the deuteron and of the pion are taken 
to be equal (r2)i/2 ^ o.657 fm and (r2)i/2 ^ 2.106 fm, re- 
spectively. The calculations were performed for the value 
of the charged pion mass M^r = 139.57018 MeV. In addi- 
tion, the calculation of the finite-size correction has been 
performed by using the latest data for the charge radii 
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(rl) = 0.452 ±0.013 fm^ [33] and (r^)i/2 = 2.1303(66) fm, 
see [34] and references therein (the result changes slightly) . 
The results of our calculations and the comparison to 
the results of Ref. [1] are given in table 1. Note that we 
have not calculated higher-order (next-to-next-to-leading) 
isospin-breaking corrections that are given in the last two 
entries of this table. The results of the calculations from 
Ref. [1] in these cases should be taken at face value. It 
can be immediately seen from the table that our calcula- 
tions completely confirm the results of Ref. [1] at next-to- 
leading order - the agreement between the two columns is 
perfect. 

After having subtracted the calculated electromagnetic 
contributions from the measured transition energy, one 
finally arrives at the strong shift, which is related to the 
nd scattering length. Since in the p-statcs the strong shift 
is proportional to and is thus tiny, the measurement of 
the quantity E^p — Eis yields directly the strong shift in 
the Is-state. In next-to-leading order in isospin breaking, 
the strong shift for the states with ? = is given by 



27rA/^n3 



1 - 



Sn{o) ~ 2 [ ^/'(n) — "0(1) -+ In a — Inn 



0(x) = r'{x)/r{x) , 



(2.4) 



where the quantity %rd denotes the threshold scattering 
amplitude in the presence of photons, which is obtained 
from the conventional amplitude by subtracting all singu- 
lar contributions at threshold (see [23-29] for more details 
and definitions). The normalization of this quantity is cho- 
sen so that in the absence of the isospin-breaking effects, 
it reduces to the ird scattering length 



%d^^T^\ 1 



and 



(2.5) 



where the ellipses stand for terms vanishing at a = and 
rUd = TTiu- These terms can be in principle evaluated in 
ChPT in a systematic manner, in analogy with more sim- 
ple cases oi ttN [35,36] and NN [37] scattering. Further, 
the quantity ^^'^ = 2d"!f'„(0)/tf'„(0) stands for the vacuum 
polarization correction to the strong shift ((5!fVj(0) denotes 
the correction of the Coulomb wave function (0) at the 
origin due to the vacuum polarization effects) . This correc- 
tion was evaluated in Ref. [31] only for the ground state. 
However, the approach used in this paper can be straight- 
forwardly generalized for the radially excited states. Fi- 
nally, it is interesting to note that the n-dependence of the 
correction term in Eq. (2.4) is universal, since short-range 
effects are the same in all atomic states. For this reason, 
even potential models (see, e.g. [32]) agree with our result 
in what concerns the difference of the correction terms in 
the states with a different n. 



Calculated corrections 
to E3p - El, [eV] 


Ref. [1] 


This work 


Point nucleus (Klein-Gordon) 


3074.69 


3074.69 


Nuclear and vr" finite size 


-0.51 


-0.52 






-0.53 [33, 34] 


Vacuum polarization a{Za) 


3.72 


3.72 


Relativistic recoil 


-0.02 


-0.02 


Higher order radiative 


0.04 




corrections 






Nuclear polarization 


0.03 





Table 1. Comparison of the electromagnetic corrections to 
the Esp — Els transition energies, calculated in Ref. [1] and in 
the present paper. The second entry for the finite-size effect 
has been obtained, using the latest experimental data on the 
charge radii. We did not address the calculation of the last two 
entries in this table. 



To summarize, in this section we have checked the va- 
lidity of the procedure which is used for the theoretical 
analysis of the pionic deuterium data at PSI. The cal- 
culated electromagnetic shift agrees very well to the one 
given in Ref. [1]. Further, we have obtained the general 
expression for the (complex) strong energy shift of the 
pionic deuterium in the next-to-leading order in isospin 
breaking, in terms of the ird threshold scattering ampli- 
tude. This relation should in principle be used to replace 
the lowest-order formula (1.1) in the data analysis. Note 
however, that the Coulomb correction which is explicitly 
displayed in Eq. (2.4) (second term in the brackets), is of 
order of 2 • 10"'^, if one replaces %rd by Eq. (2.5) and uses 
the value of the scattering length given in Eq. (1.2). Note 
also that 10"^ • • • lO'^ is an expected order of magnitude 
for the vacuum polarization contribution in Eq. (2.4), see 
Ref. [31]. Since there are no obvious reasons for having 
an anomalously large isospin-breaking correction in the 
quantity T^d either (see e.g. [1,2] and references therein), 
in the following we do not consider isospin-breaking cor- 
rections to the energy shift at all and concentrate on the 
lowest-order relation Eq. (1.1). If it turns out that the de- 
termination of the irN scattering lengths from the analysis 
of the pionic deuterium data can be performed at a few 
percent level that requires the inclusion of the isospin- 
breaking effects in Eq. (1.1), one can always go back to 
the relation Eq. (2.4). 



3 Heavy-pion effective theory for nd 
scattering 

3.1 The Lagrangian 

The findings of Ref. [20], as well as the earlier work on 
the subject (see, e.g. [14]) serve as a clear indication of 
the fact that the chiral counting is not the most suit- 
able one to be applied for the description of low-energy 
7rd scattering. Most straightforwardly, this can be visual- 
ized by comparing the contributions from the individual 
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diagrams, depicted in Fig. 2 of [20], which is reproduced 
here, in Fig. 1. The contribution from the diagrams Ib+lc 
is by two orders of magnitude smaUer than the contribu- 
tion from la, although all three diagrams emerge at the 
same chiral order. The reason for this difference is that, 
in contrary to la, the diagrams Ib+lc describe processes 
with the virtual pion emission/absorption (additional sup- 
pression at a small momenta is caused by the presence of 
the 75 -vertices in the diagrams lb and Ic). 

The above discussions lead to the conclusion that it 
will be convenient to describe the ird scattering at thresh- 
old in a framework in which the absorption and emission of 
hadrons does not appear explicitly at the level of Feynman 
diagrams, but is included in the couplings of the effective 
Lagrangian. In this manner, we arrive at the theory which 
must be in a spirit similar to HP EFT (see [19] and refer- 
ences therein). Below we dwell on some differences which 
exist between the approach used here and in Ref. [19]. 



i) The HP EFT of Ref. [19] uses the notion of the diba- 
ryon field, whereas we work in terms of the elementary 
nucleon constituents and sum up all interactions in the 
NN subsystems. After substituting the expression for 
the couplings of the Lagrangian in terms of the observ- 
ables (coefficients in the effective range expansion) , the 
"deuteron propagator" in the present paper coincides 
with the dibaryon propagator of Ref. [19] in the limit 
of vanishing effective range. We opt to work in terms 
of nucleon field in order to make the comparison with 
ChPT in the Weinberg picture more transparent. 

ii) It has been argued that the technique based on the 
introduction of the dibaryon field enables one to effec- 
tively sum up all potentially large contributions com- 
ing from the large scattering length and the effective 
range, whereas higher-order terms can be treated per- 
turbatively. The results of the present paper are ob- 
tained under an additional assumption that the effec- 
tive-range term is small, leading to some technical sim- 
plifications. This assumption is, however, not critical 
- the inclusion of the effective-range term is straight- 
forward in our approach and does not affect the con- 
clusions. 

iii) In Ref. [19] the authors have studied one particular di- 
agram that corresponds to the double-scattering con- 
tribution the the nd scattering length. In this paper, a 
systematic expansion of the quantity a^d is performed 
in the small parameter x = up to and including 
terms of order x. At this order, there are additional 
diagrams apart from the one mentioned above. 

iv) The most important question is the convergence of the 
series for the ird scattering length. We believe that 
there are (indirect) indications which testify in favor 
of the convergence. First of all, at the momenta 7 ~ 
45 MeV, the pionless effective theory gives still a rea- 
sonable description of the NN sector. As was men- 
tioned, this fact is in agreement with the observation 
made in Ref. [20] that the "modified power counting" 
in the ttc? scattering length works much better than 



Ml Ml 

\ 

a ^ b c ' 

Fig. 1. Leading-order three-body graphs which contribute to 
the vrd scattering length in the Weinberg scheme (same as in 
Ref. [20]). Shaded blobs stand for the deuteron wave function, 
and the solid and dashed lines denote nucleons and pions, re- 
spectively. Numerically, graphs b-l-c are suppressed by two or- 
ders of magnitude as compared to the graph a. 

the original chiral counting^. Yet another justification 
of the method is provided by the well-known fact that 
in the Faddeev approach, the multiple-scattering series 
for the threshold ird scattering amplitude are rapidly 
convergent, since the ttN scattering lengths are much 
smaller than the deuteron radius (see e.g. [16] and ref- 
erences therein). 

After these preliminary remarks, let us consider the 
Lagrangian of our theory. By construction, the Lagrangian 
does contains only vertices with the same number of the 
incoming and outgoing pions and nucleons. Restricting 
ourselves to the non-derivative couplings, one may easily 
write down the most general form of this Lagrangian (in 
addition, we omit below also all three-body non-derivative 
terms which contain 7r+, ttq and/or the NN-paii in the ^Sq 
state: such terms do not contribute to the Tr~d threshold 
scattering amplitude at the accuracy we are working) 

+^^(*^*--+^ + ^ + ---)^ 
+^^d+{n^-7v) + id^{'K^ X tt))?/) 

+Co(^^PaV^)UV^^PaV') + CM'^ P^i^)\i^'^ P^i^) 

+/o(^^F.V)^(V'^-P^V')^-^- + • • • , (3.1) 

where the ellipses stand for the omitted three-body terms, 
as well as for the higher-order terms in the derivative ex- 
pansion. The non-relativistic pion and nucleon fields are 
defined as tt = (771,772,773), where \/2'k± = 7ri ^ i7r2, 

7ro = 7r3 and ijj = Further, Pa and Pi denote the 

projection operators onto the ^5*0 and '^S'l states, respec- 

^ And vice versa, one may treat the HP EFT, as the system- 
atic field-theoretical realization of the counting 7 ^ <C 
47r_F^, which is heuristically implemented in the "modified 
power counting" of Ref. [20] . 
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tively 



Pn 



V8 



P, 



V8 



(3.2) 



where a and t are the Pauli matrices in the spin and 
isospin space, respectively. Note that we have not intro- 
duced an elementary deutcron field in the Lagrangian. In 
our approach, the deuteron emerges as a bound-state pole 
in the Green functions after the non-perturbative resum- 
mation of the lowest-order four-nucleon vertex. 

In the above Lagrangian, (i±,co.i,/o stand for the ef- 
fective low-energy couplings. These should be determined 
from matching of the various observables. Using dimen- 
sional rcgularization for calculating the loops enormously 
simplifies these calculations: as it is well known, all two- 
particle bubbles vanish at threshold and the results of the 
tree-level matching in the two-particle sectors remain in- 
tact. For example, the constants d± are related to the ttN 
scattering lengths through 



2n{m + M^) 



(3.3) 



and this relation remains unaffected by loop corrections. 
As concerning the constant ci, wc find it more conve- 
nient to perform the matching in the ^Si channel for the 
deuteron binding energy, and not for the scattering length 
in the np system. The difference between these two meth- 
ods shows up at higher orders. 



3.2 The deuteron 

In the two-nucleon sector of the HP EFT, there is no trace 
of pion-nucleon interactions: NN scattering is described 
completely in terms of contact four-nucleon interactions. 
The only possible loop diagrams are the s-channel bubbles 
containing the vertices with the couplings co.i. At higher 
orders, one should also include the derivative four-nucleon 
vertices. 

Consider the following connected four-point function 
in D dimensions 

(2^)^5^(pi +P2-qi- q2)GN{P;p,q) 



•4 / jD iD ,D iD 

=1 I a xid X2d yid ?/2 c 



ipiXi+ip2X2-iqiVi-iq2V2 



x(0|rV'(a;i)V(.T2)^t(y^)^t(y2)|0), 



(3.4) 



where all spin and isospin indices have been suppressed, 
and the CM and relative momenta are defined as 

P = Pi + P2 = 91 + 92 , 

P=^(Pi-P2), 9=^(91-92)- (3.5) 

The center-of-mass (CM) frame corresponds to 
P'^ = (P", 0), and _D ^ 4 in physical space. 



At the energy P" Pb(P) = Md + 'P'^ /2Md, the 
four-point function Eq. (3.4) develops a bound-state pole 
corresponding to the deuteron 



3 j* 

Gn{P]P, 9) ^ X! po'f^' ^nn + regular terms , 



-P0(P)-P0 



(3.6) 



where the deuteron wave function is defined as 



^,{p) = J d^xe*P-(0|T^(|)^(-|)|P,z), 
^J{q)^ y'dV-'9V(5,,|2^^t(|)^t(„|)|o), (3.7) 

and the sum in Eq. (3.6) runs over the polarizations of the 
deuteron spin. 

On the other hand, one may evaluate the 4-point func- 
tion Eq. (3.4) with the use of the Lagrangian Eq. (3.1) 
that amounts to the resummation of the geometrical series 
corresponding to the s-channel bubbles with four-nucleon 
vertices. Further, since at the leading order the deuteron 
is a purely '^S'l state, we may put cq = in order to get 
the deuteron pole. As the result of this resummation, one 
gets 



Gw(P;p,9) = ^ 



2P 



^(z«(pi)-p0)(«;(p2)-p^) 
ici 2Pj 



1 - ci J(P0, P) (u;(qi) - q',){w{(i2) - 92°) 
-|- terms with cq , (3-8) 
where, w{p) = m + p^/2m, and 

d^l 1 



(27r)^i {w{l) - P){w{P - 1) - P" + 
mir{l~d/2) ( ^ „n P2 ' 



(47r)3 



2m-P" + - iO , (3.9) 

4m ' 



with d~D~\. In the CM frame P = the denominator 
in Eq. (3.8) develops a pole at P° = = 2m - e. This 
gives (in d dimensions) 

f47r)^e-'-~^ 

1 - ci J(M<i, 0) = , ci = \J . (3.10) 

may (1 — c(/2) 

Finally, in the CM frame the behavior of the Green func- 
tion near the deuteron pole is given by 



Gw(P;p,g)^^ 



2P7 



iZ 



(«;(pi) - pX){w{p2) - p\) Md - P" 

2P: 

' regular terms , (3.11) 



("'(qi)-9?)(«^(q2)-92") 



8 



Ulf-G. Meifiner et al.: The pion-nucleon scattering lengths from pionic deuterium 



where the deuteron wave function rcnormaUzation con- 
stant is given by 



(47r) 



(3.12) 



mtr(2-d/2) 
3.3 Pion-deuteron scattering 

The pion-deuteron scattering amplitude can be extracted 
from the 6-point connected Green function 

i2n)^5^{pi+P2+P3-qi-q2-q3)GiP,Q;p,q;p3,q3) 



X e 



ipiXi+ip2X2+ip3X3~iqiyi-iq2y2-iq3V3 



X {Q\Ti,{x^)^{x2)n^{x3)^\y^)ij\y2)TTl{y3mc, (3.13) 

where the CM and relative momenta of nuclcon pairs are 
again given by Eq. (3.5), and P3,q3 denote the pion mo- 
menta. Near the mass shell P° Md + P^/2A/d, 
Md + Q'^ /2Md, the six-point function (3.13) develops the 
double deuteron pole. Since we are interested only in the 
threshold scattering amplitude, we may take from the be- 
ginning P = Q = and Pa = qs = 0. Then, in the vicinity 
of the pole, one has 



G{P,Q;p,q;P3,q3) X! 



Md - P° 



G, 



i^]{q) 



Md 



-\- regular terms , 



(3.14) 



with 



(27r)^<5^(P+P3-Q-g3)G., 



d^xsd^ys 



xe'P3.^-'«-'^MS,z|T^-(x3)7rl(y3)|S,j) 



(3.15) 



The residue of the quantity Gij on the pion mass shell 
yields the threshold ttc? scattering amplitude 

hm 2M, {M^ - pI){M^ - ql)G,, = i6,, %d -(3.16) 

On the other hand, in the theory with the Lagrangian 
Eq. (3.1) the 6-point Green function for vanishing 3-mo- 
menta can be given in the following form 



GiP,Q;p,q;P3,q3) 



2Pj 



(7n-pO)(m-pO)(A/,-pO) 



i(2^)^5^(P + P3-Q-93)%(^',0;P3,'?3) 



(l-ciJ(Afd,0))2 
2P, 



{m-qf){m-qO^){M^-qO) 



(3.17) 




-ca 





Fig. 2. Definition of the "truncated" Green function Rij, 
Eq. (3.17). The solid and the dashed lines denote nucleons and 
pions, respectively. 



The quantity R in Eq. (3.17) corresponds to the "trun- 
cated" Green function, and the factor (1 — ciJ{Md, 0))~^ 
emerges after the resummation of the NN bubbles in the 
■^5*1 state before the first and after the last interaction of 
the pion with one of the nucleons (see Fig. 2). Finally, 
Gi denotes the sum of all diagrams in which the virtual 
scattering of the pion on one of the nucleons occurs be- 
fore (or after) all NN interactions, or in which the first 
(or last) NN interaction happens in the ^Sq state (the 
corresponding vertex is proportional to cq). This class of 
the diagrams docs not develop a double deuteron pole, and 
contributes only to the regular part of the Green function. 
Consequently, from the comparison of Eq. (3.17) to Eqs. 
(3.14), (3.15) and (3.16) one may read of the scattering 
amplitude at threshold 



Sij%d = AfRtj{P, Q;P3, 93) 
where 



pO=QO = Md, pl=q°=M^ 



,(3.18) 



N ■ 



'2M^Z = 



Y r{i-d/2) 



(47r)f l-d/2 



2M^. (3.19) 



Hence, the prescription for calculating the threshold ttc? 
scattering amplitude is formulated as follows: in the con- 
nected 6-point Green function Eq. (3.13) omit all Feyn- 
man diagrams, where the very first or very last interaction 
occurs between the pion and nucleon, or between NN- 
pair in the ^Sq state. In the remaining diagrams, resum 
all initial- and final-state NN bubbles and write the fi- 
nal result in a form of Eq. (3.17); read off the quantity 
Rij{P,Q;P3,q3); perform the mass-shell limit, let all 3- 
momenta vanish, multiply by the normalization factor Af, 
given by Eq. (3.19) and get the threshold scattering am- 
plitude %rd- 



3.4 Leading order 

At the lowest order in the expansion parameter x = j /M^^, 
there is a single contribution to the quantity Rij defined 
by Eq. (3.17), which is depicted in Fig. 3. At threshold, 
this contribution equals to 



R 



(0) 



icjd+TriP^P,] 



thr 



d^l 



1 



(27r)^i {w{l)~P)^ 



1 



2 J 3-1 



w{l) - Md + /o 



An 



(3.20) 
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P-l 



Fig. 3. Leading-order contribution to the vrd scattering length. 



Substituting this result into Eq. (3.18), using Eqs. (2.5), 
(3.3) and (3.19) and expanding = 2m + 0{e), at the 
leading order we finally obtain 



1+AU/m 
1 + AU/2Tn 



2a+ + 0{x) , 



(3.21) 



which of course coincides with the well-known result. Here 
we only wish to note that our result is valid at all orders 
in the chiral expansion for the scattering length a+. On 
the other hand, if one works in the Weinberg scheme, one 
has to identify the contributions to the quantity a+ order 
by order in the chiral expansion [20]. 



3.5 Next-to-leading order 

Since the quantity a+ turns out to be very small, the cor- 
rections exceed in magnitude the leading-order result. In 
the HP EFT, four diagrams depicted in Fig. 4 contribute 
at next-to-leading order. At threshold, we get 



-Rij'' =5ijc\{Ra + Rb + Rc + Rd) , 



Ra^2{d\-2<f_] 



1 



(27r)^^ {2n)Di w{l) - 1° 
1 



iwi\) - Aid + mi" - fc" - Md + (1 - k)2/2M,) 
1 



(w(k) + fcO)(u'(k) - Md - fcO) ' 



^^^2(^ + 2^-)/ (2^(2^ 



1 



(27r)^^ {2n)Di w{l) - P 
1 



(w(l) - Md + m{w{l - k) - Md - M^ + P- k") 
1 



Rc = M 



d^l d^k d'^q 1 



+ 7 (27r)«z (27r)^z (27r)^i w{\) - P 
1 

{w{\) - Md + ?")(«;(! - q) - - M^ + - qO) 

Cl 1 

1 - ciJ{Md + M^ + qO, -q) AU + 9° + qy2AU 



w{k - q) - Aid - AU + k° - q° 



(u;(k) - fcO)(w(k) - Md + k°) 



1 



Rd = ^ fa 



(3.22) 



Note that performing (formally) the limit Ali^/m in 
the quantity one gets l/{f-k"-Md+{\-'kf/2M^) 
2Mt^/ (1— k)2. In this limit, it is possible to relate the quan- 
tity Ra to the average of the operator l/q^ between the 
deuteron wave functions in the potential that corresponds 
to pointlikc interaction ci{ip'^ Pitp)'^ {tp'^ Pi^). In the same 
normalization as in Rcf. [20] one obtains 



-Ra ^ 2{d 



2dt 



2AU 1 
Z {2tt) 



(3.23) 



w.f . 



and the standard expression for the double-scattering con- 
tribution in the limit M.^/m — *■ (see e.g. [20]) is repro- 
duced. 

The counting of the above diagrams proceeds as fol- 
lows. According to Eq. (3.10), the coupling ci counts like 
cc^^, and the couplings d± count like x" . Further, after in- 
tegrating over the time-like components l'^\k" ,q" , one may 
rescale I ^ 7I, k ^ 7k, q — *■ 7q, with 7 = 0{x). Each 
propagator of a pion or a nucleon counts as 7^^ ^ a;-2 
and the "virtual deuteron propagator" ci/(l — ciJ) counts 
as Cl ~ x^^. With these counting rules, it is straightfor- 
ward to ensure that Ra ^ Rb ^ Rc = 0{x") (modulo 
logarithms). Furthermore, since the constant /o cancels 
the ultraviolet divergences in the diagrams 4a,b,c, it must 
count at the same order in x. This fixes /o = 0{x~'^). 
Note also that the contributions proportional to the cou- 
pling constant Cq (nmr'^ intermediate states) have been 
dropped from Rc altogether at this order. This is related 
to the orthogonality of the projectors Pi and Pa given in 
Eq. (3.2). 

The fact that the quantity Rc in Eq. (3.22) is propor- 
tional to simplifies the calculations considerably. Since 
a+ is very small, we shall systematically neglect in 
all expressions, and thus assume Rc = 0. Evaluating the 
remaining integrals in dimensional regularization and car- 
rying out the renormalization in a standard manner, we 
finally arrive at the following expression for the ird scat- 
tering length at the next-to-leading order 



(1) _ m(l + M^/mf 



7r(l + M^/2m) 



X a_ 



{Xa - Xb) In 



rae 



+ (Ja - Jb) 



0{x^ 



Alt 



A7r^{l + AU/2m) 



(3.24) 



where foip-) denotes the renormalized coupling constant 



/o 



647r3(l-FM^/TO)2 2 



eAB 



a_ {Xa - Xb)\ + /o(m) : 



10 



Ulf-G. Meifiner et al.: The pion-nucleon scattering lengths from pionic deuterium 





c 

d 

Fig. 4. Next-to-leading order contributions to the nd scatter- 
ing length, see Eq. (3.22): (a) the standard double-scattering 
contribution, (b) rescattering on a single nucleon, (c) pion 
rescattering on the infinite chain of the nucleon bubbles ( "vir- 
tual deuteron" , denoted by a thick dot-dashed line) , (d) coun- 
terterm contribution. 

and ^ denotes the scale of dimensional regularization. Fur- 
ther, Xa,bj Ja,b denote the integrals over Feynman param- 
eters which depend on the dimensionless variable MT^/m 
and emerge from Ra,b- These integrals are evaluated in 
appendix B. Here, we only give their approximate values 

xa = -1.2569, Ja = -1.1334, , . 

Xb = 0.5098, Jb = -0.3731 . ^ > 

At present, the numerical value of the countcrtcrm fQ^fJ.) 
is not known. For this reason, one has to include this un- 
known quantity completely in the theoretical error and to 
estimate the uncertainty that emerges already at next-to- 
leading order. Most easily, this can be done by using the 
renormalization group equation for the scale dependence 
of /o which at this order reads 

d fn N 87r(l + M-^/mf ^ 

^dil-^^^^'^ 7m^ ^- ^ ' 

We use the following procedure to estimate the uncer- 
tainty. Since the hard scale of the theory is of order of 
M.^, we may set /o(m) = at a some scale ~ M^r and 
study the /^-dependence of the plot in the (a+, a_)-plane, 
which emerges from Eq. (3.24) at a given (experimental) 
value of a^d. This plot is given in Fig. 5. Varying /i in a 
"reasonable" range, we may thus visualize the error that 
is caused by the scale dependence. Here, we wish to note 
that the scale dependence is of course not the only pos- 
sible source of the theoretical uncertainty in general. In 
order to have a reliable estimate of the error (in the case 
of the weak scale dependence) one has, in addition, to use 
dimensional arguments to estimate the size of the LECs. 




a. in units of (M^^) 



Fig. 5. The constraints on the S-wave nN scattering lengths 
a+ and a_ from the pionic deuterium energy shift for different 
values of the scale parameter fj,. The central line is given for 
/i = 146 MeV. The upper and lower bounds correspond to 
jj, — 250 MeV and = 100 MeV, respectively. In addition, 
we display the constraints from the most recent measurement 
of the pionic hydrogen shift and width [38]. The upper limit 
for the width is measured in the 4p — Is transition in pionic 
hydrogen. 

However, in the case of a strong scale dependence, as in 
the example considered here, additional arguments are not 
needed. 

The results which are displayed in Fig. 5 are in a qual- 
itative agreement with the findings of Refs. [18, 19]. Note 
that these results arc obtained at the ncxt-to-leading or- 
der in HP EFT. It is unlikely that taking into account 
higher-order terms will reduce the uncertainty due to the 
scale dependence. On the contrary, from the comparison 
with e.g. Eq. (18) of Ref. [19], one may conclude that 
numerically the most important corrections due to the ef- 
fective range at this order amount to the multiplication 
of the loops by the deuteron wave-function renormaliza- 
tion factor Zd ~ {I — "yrd)~^ ~ 1.7, where = 1.765 fm 
is the parameter related to the effective range in the ^Si 
channel. This effect leads to further amplification of the 
ambiguity related to the scale dependence. 

The interpretation of the results displayed in Fig. 5 is 
unambiguous. The experimental data together with the 
above theoretical analysis constrain the 5'-wave ttN scat- 
tering lengths within the band which - for large values 
of /i - also intersects with the common area of two other 
bands, obtained from the data on the pionic hydrogen en- 
ergy shift and width. However, without having in advance 
estimated the numerical value of /o, one can not use the 
measured value of aT^d for a precise determination of the 
scattering length a+, that was an original motivation for 
studying the pion-deuteron system. Moreover, since the 
constant fo parameterizes the short-range physics, it can 
be only fixed either by using other experimental data (dif- 
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ferent from the Trd scattering process considered here) or 
through the lattice simulations. 

We wish to add that, in order to be consistent with 
power-counting, the magnitude of the LEG /q (/i) is para- 
metrically enhanced by a factor Xjx^ ~ 10 as compared to 
the dimensional estimate /g (^f) = const/M|. Numerically, 
substituting the dimensional estimate into Eq. (3.22), one 
gets the uncertainty ~ 3% in the quantity a-^d^ which 
by far underestimates the actual uncertainty displayed in 
Fig. 5. 

3.6 Imaginary part 

In the presence of the absorptive channels, the coupling 
constant /o is not real (note that in order to simplify the 
notations, we have up to here always omitted the symbol 
"Re" in this coupling constant, as well as in a-j^d and other 
quantities, if this docs not lead to the confusion). The 
imaginary part of the -nd scattering length is given by (cf. 
with Eq. (3.22)) 

Ima.d= ^^glm/o, (3.28) 

where ^\ = ^'^ /2ti denotes the square of the dcuteron 
wave function at the origin. Further, the imaginary part 
of the constant /o is directly related to the inelastic chan- 
nels, which are "shielded" when constructing the HP EFT. 
Since we have neglected all electromagnetic effects from 
the beginning, the bulk contribution to Im /o is provided 
by the two-neutron intermediate state, whereas the con- 
tribution from the ^nn intermediate state is omitted. To 
take this fact into account, wc add the superscript "str" to 
the pertinent scattering length. To get the imaginary part, 
which also includes the effect of the ^nn state, one writes 
Im a,rd = Im a^'^ (1 + 1 /-R7) i where the experimental value 
of the quantity = 2.83, which stands for the ratio of 
the cross-sections of the ■K~d transition into the nn and 
7nn final states (the Panofsky ratio) and is treated at face 
value. At the lowest order we get 

X {2^f5\pi + P2 + P3 - 91 - 92), (3.29) 

where T(pnn^ nn) stands for the transition ma- 
trix element for the process p{pi) + n[p2) -\- t^~{p3) — > 
n(fci) -I- n{k2), and the above integral has to be evaluated 
at threshold Pi = P2 = P3 ^ 0. After straightforward 
calculation we get 

Im/f = ^|T(pn7r-^nn)p, 



|ki| = |k2| =p* = ^mM^ + --- , (3.30) 

where the ellipses stand for the relativistic corrections. 
Now we note that the above result is compatible with the 



well-known relation in terms of the n d ^ nn inelastic 
cross-section at threshold 

Ima^d = T- lip IPsI cr(7r^d nn) , 
47r IpaHo 

a{TT~d nn) ~ ^'^^ |T(7r~o? nn)p , (3.31) 

if the transition amplitudes for the processes ir^d — > nn 
and pmr~ nn are related by 

Tin^d^ nn) =<PoT{pmr- ^ nn) . (3.32) 

The equation (3.32) is nothing but the leading-order Deser 
formula, which is obtained with the assumption that the 
deutcron radius is much larger than the distances rele- 
vant for the interactions in the ttNN — NN system. The 
corrections to this formula would then emerge at 0{x) rel- 
ative to the leading-order result. Note also that Eq. (3.32) 
agrees with the result of [13], obtained within the poten- 
tial scattering theory, that serves as a good check for the 
vahdity of HP EFT. 



4 Matching to ChPT and numerical analysis 

4.1 Threshold amplitude in the HP EFT 

As it was demonstrated in the previous section, within the 
HP EFT it is not possible to get an accurate description 
of the nd scattering length in terms of the S-wavc ttN 
scattering length only. The relation between these quan- 
tities contains a large unknown short-range contribution 
(three-body force in the language of the potential scat- 
tering theory), which is parameterized through the LEG 

In this section we shall address the issue whether it is 
possible to achieve an increased accuracy if one treats the 
same problem within the Weinberg approach [14,20]. In 
the HP EFT (which is the effective theory of ChPT in the 
Weinberg picture for the momenta p <C Mjr), the LEGs 
(including /o) receive contributions from two different mo- 
mentum regions: M^r < p < A and p > A (here A ~ 47rF^ 
denotes the cutoff mass used in the Weinberg formulation, 
which is of order of the hard scale in ChPT) . In the Wein- 
berg framework, one has "resolved" the momenta at the 
scales < p < A. The unknown dynamics at the mo- 
menta p > A is parameterized by new LECs which are 
now defined at the scale A instead of M^. On dimensional 
grounds, the natural size of these new LEGs must be much 
smaller than the size of old LECs, since M^r ^ A. Stated 
differently, if we start from the LECs in the theory at a 
scale A and calculate the LECs of the HP EFT in the limit 
of a large M^, we must see that the LECs of the HP EFT 
must be enhanced by the pion loops where the loop mo- 
mentum runs within < p < A. It is natural to assume 
that this momentum region contributes the bulk of the 
total magnitude of the LEG in question. Could one then 
separate the large but potentially calculable pion exchange 
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contribution to /o from unknown short-range contribution 
at a scale A? 

A natural choice of the 5-matrix element in the scat- 
tering sector, which can be used to determine the con- 
stant /o through the matching procedure to the Weinberg 
framework, is that of the elastic process tt~ {pn)3g^ — > 
Tr~ {pn)3g^. The matching condition has the form 

3 

Tw{piP2P3;qiq2q3) = l[m{p,)]'/^m{(i,)]'/^ 

1=1 

xT{piP2P3;qiq2q3) , (4.1) 

where T and Tw denote pertinent scattering matrix ele- 
ments in the HP EFT and in ChPT, respectively. Further, 
E,{1) = + P, i = 1, 2 and £3(1) = ^ + P are the 
relativistic energies of a nucleon and a pion. Note that the 
kinematical factor in Eq. (4.1) is introduced in order to 
take into account the different normalization of the one- 
particle states in HP EFT and in the Weinberg approach. 

We found it convenient to perform matching for the 
one-particle irreducible (IPI) matrix elements separately 
in two- and three-particle subsystems (with respect to 
the non-relativistic pion and nucleon propagators). In the 
three-particle sector, the scattering amplitude turns out to 
be singular at threshold: one has to choose the particular 
kinematics in order to approach the the zero-momentum 
limit. A possible choice of the external momenta is given 
by 

p = q = 0, P3 - -P, qs - -Q, P - -Q ^ 0, (4.2) 

where the definition of the CM and the relative momenta 
is given in Eq. (3.5). 

Below, we schematically describe the matching for the 
above scattering matrix element. Explicit expressions will 
not be used in the discussion. For our purposes, we only 
need to demonstrate that such a matching can be per- 
formed in principle. 

We start with the evaluation of the scattering ampli- 
tude T that enters the right-hand side of Eq. (4.1). This 
quantity receives many contributions (some of them are 
depicted in Fig. 6). Out of these contributions, the dia- 
grams 6a and 6cl are one-particle reducible and will be 
excluded from matching. Some of the contributions are 
singular at threshold. For example, the diagram in Fig. 6e2 
contains the factor |P|~^, and the diagrams in Fig. 6e3,e4 
are logarithmically singular at small momenta. The three- 
particle threshold scattering amplitude A{fJ.), by defini- 
tion, is obtained from the IPI matrix element, by subtract- 
ing first all singular contributions at |P| — > 0. Namely, in 
the vicinity of threshold, the sum of the IPX diagrams 
shown in Fig. 6 has the form 

T(iPi) = ^ To In ^ + Ait,) + 0(|P|) , (4.3) 
|P| M 

where the arbitrary scale that enters the logarithm is set 
equal to the scale of dimensional regularization /i, for sim- 
plicity. Further, the quantity A{^.) = /o (m) + A{^.) con- 
tains additive contribution from the LEC foip-), which 




a b 



+ /' + \ / + 



cl c2 c3 




el e2 63 e4 



Fig. 6. Lowest-order diagrams contributing to the threshold 
scattering amplitude for the process tt~ (pn)3g^ — > n~ (pnjsg^. 
The solid, dashed and dot-dashed lines correspond to nucle- 
ons, to pions and to the deuteron (an infinite sum of nucleon 
bubbles), respectively. All diagrams except (a) and (cl) are 
one-particle irreducible. 



emerges in the diagram Fig. 6b. Note that the corrections 
to this contribution due to the initi al- an d finite-state in- 
teractions vanish in HP EFT. Here, .4(/i) stands for a sum 
of all diagrams that do not contain /q(/x), see Fig. 6. At 
the order we are working, this quantity depends only on 
the couplings in the two-particle sector d±, ci, as well as 
the masses M^, m. 

The power counting in the expansion parameter x is 
organized as follows. One groups the diagrams shown in 
Fig. 6, according to the number (djv) of virtual A^A^ inter- 
actions - "squeezing" all deuteron propagators to a single 
point. For example, the IPI diagrams (c2),(c3) correspond 
to dN = 0, (dl), (d2) - to dN = 1, (el), (e2), (e3), (e4) - 
to dN = 2, (fl) - to d^v = 3, (gl) - to = 4, etc. 

For illustration, let us first consider the IPI diagrams 
with d^ = 0. The contribution to the threshold amplitude, 
which is obtained after subtracting all singular pieces from 
these diagrams, is proportional to powers of the ttA^ cou- 
pling constants d± and depends on the masses M^r, m and 
the scale fi. Neither of these quantities scale with x and 
thus this contribution emerges at 0(1). 

Next, we consider the diagrams with djv = 1. The dia- 
gram (dl) is proportional to cid^. = 0(x~^). In addition, 
there is an intrinsic scale x present in the deuteron prop- 
agator (1 — Cl J)~^, since ci = 0{x~^). In order to get rid 
of this scale, one has to first perform the contour integra- 
tion over the time-component of the loop momentum k'^ 
and afterwards rescale the loop three-momentum k — > zk. 
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as well as all external three-momenta. Further, after sub- 
tracting all singular pieces at threshold, one may put the 
external momenta to zero. The counting of the powers of 
X proceeds as follows. The integration measure d'^k yields 
the factor . Since the integration over reduces the 
total number of elementary (pion and nucleon) propaga- 
tors from three to two and since each elementary propa- 
gator counts like x~'^ , this leads to the factor after 
rescaling the momenta. Finally, taking into account the 
fact that after rescaling the deuteron propagator counts 
like 0(1) and putting together all factors, we arrive at the 
conclusion that the diagram (dl) counts like x~^ (modulo 
logarithms). By using the same method, it is easy to show 
that the diagram (d2) contributes at 0{x~^). 

Applying the same argument to all diagrams in Fig. 6, 
we finally come to the conclusion that the leading order 
scaling of the loop diagrams is that at 0(x~^). Only the 
diagrams (dl), (el), (e3), (e4), (fl) contribute at the lead- 
ing order in x. All other diagrams, as well as the diagrams 
which can be obtained from the diagrams displayed in 
Fig. 6 by attaching more pion and/or nucleon loops, con- 
tribute at a higher order in x. 

Finally one has to consider the LEG /o(/^) which, as 
we know, scales at 0(x~^) at leading order. Generally, one 
can write 



r^i^i) = F^2X-^ + F^^x- 



Fo + Fix - 



(4.4) 



where only the leading coefficient F^2 contributes in the 
pion-deuteron scattering length at 0{x). 

To summarize, the right-hand-side of Eq. (4.1) can be 
written as a sum of the terms that scale as with n = 
— 2,-1,0, The leading-order contribution at 0{x~'^) 
emerges from the diagrams (b), (dl), (el), (e3), (e4) and 
(fl). Except the LEG /o that enters from diagram (b), all 
other diagrams are expressed in terms of the parameters 
determined in the two-particle sector. 



4.2 Weinberg- Tomozawa term 

Next, we consider the evaluation of the left-hand side of 
Eq. (4.1). The lowest-order contribution in GhPT in the 
Weinberg framework emerges from the diagram depicted 
in Fig. 7. This diagram describes the double scattering of 
the pion on the nucleons, with the vertices obtained from 
Weinberg- Tomozawa Lagrangian 



C = 



{(7r+9^7r0-^0a,,^+)p7^r 



-I- (tt^S^jTt" - 7r''(9p7r^)n7^p} 

+ ^ (7r"9,,7r+ - 7r+a^7r-){p7> - n7''?i} , (4.5) 

where one has used the Gondon-Shortley convention for 
the component fields and F stands for the pion decay 
constant in the chiral limit. At the order of accuracy we 



Fig. 7. Lowest-order contribution to the threshold amplitude 
of the process n~pn n^pn in ChPT. The solid and dashed 
lines denote nucleons and pions, respectively. 



are working, we may take F ^ F^, = 92.4 MeV. The 
scattering amplitude for the process p{qi)n{q2)T^~ {q^) 
p{pi)n{p2)TT~ {ps) at second order in perturbation theory 
is given by 



w 



1_ u{pi) 4-iu{qi) u{p2) ^3^(92 

4^4 _ _ _ ^3)2 

1 u{Pi) ^3'^{Qi) u(p2) 43u{q2) 
4F4 Mi^{pi^qi+p^Y 

1 u{Pl) ^3u{q2) u(p2) 4zu{qi) 

2F4 Af2-(pi-q2+p3)2 



(4.6) 



The low-energy reduction of the relativistic amplitude 
Eq. (4.6) yields both the reducible and irreducible parts. 
In order to single out the IPI piece, one has to split the rel- 
ativistic pion propagator into the positive- and negative- 
energy components and expand the result for small three- 
momenta. For example, 

1 



Ml - (pi - 91 - <Z3)2 
1 



AU - g? - <?§ + + (qi + qa - Pi) 

1 

4MI^'" ' 



(4.7) 



where the first (the second) term correspond to the one- 
particle reducible (irreducible) pieces. Evaluating the IPI 
part at threshold, we get 



(IPI) 



w 



2F^ 



(X1X1X2X2 + X1X2X2X1): (4-8) 



where the Pauli spinors are defined through \/2m ^ "q* ^ ~ 

lim^^^o "(9i) and -\/2m (xliO) = linip--fO u{pi). Perform- 
ing a Fierz transformation, we get 



X 



(4.9) 



Finally, dividing the amplitude Tw by a kinematical factor 
{2m,)'^2AlTr that emerges from Eq. (4.1) at threshold, one 
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may read off the value of the LEG /o, which is obtained 
through the matching to ChPT at the leading order in 
chiral expansion 



/o = 



1 



(4.10) 



Substituting now this value into Eq. (3.22), one ends up 
with the tiny contribution —0.0005 (cf. with the cur- 
rent experimental value OTrd = 0.0261 M~^). The small 
magnitude for this correction is in agreement with find- 
ings of Ref. [14]. The reason why the leading-order chiral 
contribution is so small is simple: it emerges only at the 
NNLO in the x counting and is contained in the term Fq of 
Eq. (4.4). To the contrary, some of the sub- leading contri- 
butions in ChPT, which are contained in the LO and NLO 
coefficients F_2, -F-i, get enhanced by inverse powers of 
X and are numerically much larger that the Weinberg- 
Tomozawa term. The above example clearly shows that 
the power-counting in the HP EFT and in the Weinberg 
picture are not correlated. In accordance with the findings 
of Ref. [20] , the counting which is based on the expansion 
parameter x, better reflects the numerical size of the con- 
tributions emerging at different orders. 



4.3 Initial- and final-state interactions 

From the discussion in the previous section one concludes 
that one needs to identify a sub-class of the diagrams 
in the Weinberg approach which - after the matching - 
contributes at 0{x^^) to the LEG fa. Since here we are 
only interested in establishing accuracy limits in the cal- 
culations of the pion-deuteron scattering length, in the 
matching condition we may neglect loop diagrams shown 
in Fig. 6 and their counterparts in ChPT. Since these di- 
agrams depend only on the parameters that can be deter- 
mined in the two-particle sector, including these diagrams 
will only shift the central value of /o determined through 
matching, without significantly affecting the error bars. 

The contributions which one needs to retain in the 
left-hand side of the matching condition Eq. (4.1), con- 
tain the genuine 3-particle ChPT LEC(s) /q, which are 
the counterpart (s) of the LEC /o m HP EFT. Note that 
due to chiral symmetry, the LECs in the Weinberg frame- 
work appear first at O(p^) - the corresponding Lagrangian 
should contain either two derivatives of the pion field or 
the quark mass matrix. Further, these LECs do not come 
alone: in the spirit of the Weinberg approach, one has in 
addition to consider summing up the diagrams that corre- 
spond to the initial- and final- state interactions of pions 
and nucleons. Due to the use of the cutoff in the Wein- 
berg approach, this effect does not vanish at threshold, in 
difference to HP EFT. 

In the following, we shall first provide a very crude 
estimate of the initial- and final-state interaction effect. 
We assume that the bulk of the effect in the 3-particle 
system comes from the NN pair, whereas the "hopping" 
of the pion on the nucleon lines amounts to a small cor- 
rection. We further approximate the one-pion exchange in 



-OWO- 
R *'o R 



+... 



Fig. 8. Including the initial- and final-state interactions in the 
Weinberg picture 



the Weinberg picture by a local NN potential with a cou- 
pling denoted by c^, assuming a cutoff between the pion 
mass and the hard scale A in the loops, wherever this 
vertex is present. With these assumptions, the matching 
condition takes the form that is schematically shown in 
Fig. 8 (in this figure, the ellipses stand for the diagrams 
which do not contain 3-particle LECs). One may introduce 
the "amplification factor" which is obtained by summing 
up all bubbles in the ingoing and outgoing NN lines 



^^(IpI) 



1 



{l-c[iA)JAiOW 



(27r)3 P - p2 - jO 



mA 



where we have made explicit the yl-dependence of the cou- 
pling constant in the Weinberg scheme. 

Imposing now the condition that the denominator has 
a pole at the deuteron binding energy |p| = ij, one can de- 
termine the coupling constant c'^ {A) and the amplification 
factor 



7717 

47r 



and the matching gives 

/O (X i? 2/o + • 



iA^ 



r^x 



2H 

Jo 



(4.12) 



(4.13) 



where does not scale with x. Hence, the desired scale 
behavior in the constant /o emerges from the sub-class of 
the diagrams in the Weinberg picture which describe NN 
rescattering in the initial and final states in all orders. 

The numerical value of the amplification factor 
turns out to be very large: 7?^ ~ 4 for yl = M^, 7?^ ~ 
50 for A = 500 MeV and i?^ ~ 200 for yl = 1 GeV. Per- 
forming now a standard dimensional analysis of the LECs 
in ChPT (see e.g. [39]), we get the estimate of the uncer- 
tainty in the quantity /o 



1 



2AL 



1 



2F^M, V 4^F, 



(4.14) 



Substituting this result into Eq. (3.22), one finally gets 
Z\a„rf/a^d = [0.5%; 6%; 23%] 



for A = [M^; 500 MeV; 1 GeV] 



(4.15) 
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Now we present the exact numerical results in the Wein- 
berg picture which confirm the validity of our crude esti- 
mate. Note that, if the calculations are done in the Wein- 
berg approach ab initio^ the amplification factor is con- 
tained in the quantity |f?o(^)P, where ^o(^) denotes the 
value of the wave function at the origin. The counterpart 
of this quantity in HP EFT is given by (pQ = {j^/2Tr)i = 

0.0445 fm 2 , For comparison, the NNLO wave functions 
in the Weinberg approach are given by 
^o(^) = [0.487; 0.434] fm"^ for A = [450; 650] MeV [40]. 
Hence, one gets the following estimate for the amplifica- 
tion factor i?^ = <I>o{A)^/<pQ ~ 100, which in turn cor- 
responds to the reasonable value of the cutoff mass A ~ 
720 MeV in our order-of-magnitude estimate in Eq. (4.15) 
and to AaTrd/ciTTd — 12%, which qualitatively agrees with 
the large uncertainty in Fig. 5. Note that the NN interac- 
tions are now studied up to and including N^LO [41]. In 
our estimates we however use the NNLO wave functions, 
in order to consistently compare with earlier work on the 
subject [20]. 

The LEG /o that enters the expression of the pion- 
deuteron scattering length, can be estimated by using res- 
onance saturation. The details of this procedure can be 
found in appendix C. It turns out that using the set of pa- 
rameters determined in Ref. [42] leads to the result which 
is consistent with the above dimensional estimate. How- 
ever, due to the large uncertainty in the values of these 
parameters, the resonance saturation hypothesis can not 
be used - at present stage - to get a better accuracy in 
the calculation of aT^d- To this end, one has to improve 
the quality of the fit of the parameters of the resonance 
Lagrangian to the data. 

How can this large uncertainty be reconciled with the 
mild cutoff dependence of Ref. [20]? The answer to this 
question is the following. Assuming for instance that the 
cutoff dependence of the generic LEC(s) /q in ChPT at 
the lowest order is logarithmic, we may split these LEC(s) 
into a term which depends of the cutoff mass A and a term 
which does not 

/^^A + i?lnA, (4.16) 

where A* stands for a some characteristic hard scale (typ- 
ically ~ 1 GeV). The mild cutoff dependence is equivalent 
to the statement that the constant B is small. One may 
argue that the smallness of B is related to the dominance 
of the one-pion exchange which at the large distances has 
a softer behavior than the contact 4-fcrmion vertex which 
is used in HP EFT to describe the deuteron. However, the 
constant A is not covered by above argument. Using the 
dimensional arguments, wc see that the bulk of the un- 
certainty ~ 12% should come from the scale-independent 
constant A. 



5 Comparison to the existing approaches 

The issue of the pion-deuteron scattering has obtained 
an extensive coverage in the literature during the last 



decades. A highly incomplete list of references is given 
in [5-13], see also references therein. Note that the frame- 
work used in HP EFT is in fact very close to that of the 
potential scattering theory. Consequently, a rapid conver- 
gence of the series for CTrd in HP EFT would indicate also 
on the applicability of the potential picture in the ird scat- 
tering at threshold. The main question related to choice 
of the potential remains however open. As was demon- 
strated in Ref. [43] , the conventional quantum- mechanical 
potentials can be interpreted as a mere rcgularization of 
the non-relativistic effective field theories. From this point 
of view, the scale dependence which was discovered in this 
article is equivalent to the off-shell effects in the two-body 
potentials which must be canceled by the corresponding 
three-body force (analog of the LEG fa). It is clear that 
fitting two-body potentials to the scattering data can not 
completely eliminate this inherent off-shell (scale) depen- 
dence. Another aspect of the problem concerns the ab- 
sence of the relation to QGD and to GhPT. One may con- 
clude that the existing potential models which are used 
to extract ttN scattering lengths from the measured nd 
scattering length can not in principle provide enough ac- 
curacy needed for the test of the predictions of QGD at 
low energy in the ttN sector. 

Next, we shall consider the "hybrid" approach in nu- 
clear physics, which is based on the calculation in GhPT 
of a certain set of Feynman diagrams, corresponding to 
the "irreducible transition kernel" and finally sandwich- 
ing the result by "realistic" wave functions which are cal- 
culated, using Paris, Bonn, Argonne, . . . potentials (see 
e.g. [12, 14, 16]). As we already mentioned above, this pro- 
cedure can be justified, if and only if the long-range effects 
(e.g. the one-pion exchange) dominate the transition oper- 
ator. However, as we have seen, the situation in describing 
nd scattering length is just the opposite: this quantity re- 
ceives a large short-distance contribution from the LEG 
/o, which should be there in order to provide a scale- 
independence of the final result. In the final expression 
for the decay width, this constant is multiplied by the 
wave function of the deuteron at the origin squared. Gon- 
sequently, for the consistency of the hybrid approach, one 
must be sure that the value of the "realistic" wave function 
at the origin is a good approximation of the same quan- 
tity, obtained in the effective theory which was previously 
used to calculate the transition kernel. 

A fully consistent approach to the pion-dcuteron scat- 
tering problem is provided by the effective field theories, 
in which the wave function of the deuteron is evaluated 
within the same setting as the diagrams, describing the 
irreducible kernel. These are e.g. the calculations carried 
out within the framework with perturbative pions [18], HP 
EFT with elementary deuteron field [19], or in the Wein- 
berg scheme [17,20]. The large scale dependence in the 
LEGs analogous to our /o, was first reported in Ref. [18]. 
The Feynman diagram which leads to such a large de- 
pendence is the counterpart of our Fig. 4a. One should 
however note that treating the pions relativistically and 
using the chiral Lagrangians unnecessarily complicates 
the simple physical picture, since the question about the 
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convergence of ChPT expansion naturally arises. A more 
straightforward approach is provided by HP EFT [19]. The 
authors of this article also find the logarithmic enhance- 
ment of the diagram 4a, as well as the large scale depen- 
dence which emerges due to this diagram. Note however, 
that in this paper not all terms at 0{x) have taken into 
account: The diagrams shown in Fig. 4b. c arc omitted - 
although, as we have seen, the inclusion of these diagrams 
does not change the qualitative conclusions. In the present 
paper, we have also critically re-examined the conjecture 
made in Ref. [19], concerning the possibility of calculation 
of the LEG /o with an improved accuracy in the Weinberg 
picture - this, as was shown, is not possible. 

Finally, we briefly comment on the papers [17,20], 
which provide the systematic treatment of the nd scat- 
tering problem within the Weinberg framework. The re- 
sults of these papers are now very easy to understand and 
the relation to the HP EFT becomes crystal clear. Thus, 
the puzzle concerning the (seeming) differences between 
the HP EFT and the Weinberg approach has finally been 
resolved. 



6 Conclusions 

i) In this paper, we systematically investigate pionic deu- 
terium within the framework of effective field theories. 
The whole treatment naturally falls into several steps. 
At the first step, we discuss the extraction of the ird 
scattering length from the 3p — Is transition energy 
and width in the pionic deuterium. The next-to-leading 
order result has been obtained for the level energies, 
which can be used for an accurate determination of the 
ird threshold scattering amplitude from experimental 
measurements. Since the isospin-breaking corrections 
in this amplitude are not expected to be relevant, given 
the relatively large theoretical uncertainty in connect- 
ing with the irN scattering lengths, these correc- 
tions have been neglected for the time being. 

ii) The main focus in the present paper is on investigat- 
ing the possibility to relate aj^d to the ttN scattering 
lengths and on the analysis of the systematic theoret- 
ical error in such a procedure. We give a consistent 
treatment of the problem within the framework of HP 
EFT, where the expansion parameter is given by the 
quantity x = ^/M.^ ~ 1/3, where 7 ~ 45 MeV is the 
characteristic bound-state momentum in the deuteron. 
In this paper, we have evaluated contributions to the 
quantity a-^d, up to and including 0(x) that corre- 
sponds to the next-to-leading order in HP EFT. 

iii) At next-to-leading order, the ird scattering length re- 
ceives a contribution from the (unknown) LEG which 
we denote as /q. We have used the scale-dependence 
of this contribution for estimating the theoretical er- 
ror in our calculations. This scale dependence is shown 
in Fig. 5. The scale must be chosen to be of order 
of the pion mass, but is otherwise arbitrary. In our 
opinion, the range 100 MeV < ^ < 250 MeV can be 
roughly considered as a "natural" choice of this scale. 
It should be also taken into account that the plot in 



this figure, which corresponds to the recent experimen- 
tal measurements of the pionic hydrogen decay width, 
still does not include the isospin-breaking corrections 
from ChPT [28, 44]. As we see from Fig. 5, the theoret- 
ical uncertainty due to the unknown LEG /o is rather 
large. 

iv) In this paper, we have investigated in detail the dif- 
ferences between the Weinberg approach and the HP 
EFT. It was shown that, despite the very mild cutoff 
dependence in the Weinberg approach, the uncertainty 
due to the unknown LECs is significant and is of the 
same order of magnitude as in the HP EFT. The rea- 
son for this is that the large initial- and final-state 
NN interactions lead to the amplification of the ini- 
tially small LEG contribution. Taking into account this 
amplification, the theoretical predictions within both 
approaches are essentially the same. 

v) Our main conclusion, concerning the accuracy limits in 
the extraction of the nN scattering lengths from the 
pion-deuteron data, can be formulated as follows: by 
far the largest source of uncertainty is the low-energy 
constant /o, which is the genuine short-distance three- 
body contribution and should be either determined 
by other experiments or should be obtained by lat- 
tice simulations. In particular, one might attempt to 
get at least the order-of-magnitude estimate from the 
process NN NNtttt or from the pion-nucleus opti- 
cal potential, in case of non-equal proton and neutron 
densities, where /o should be present apart from the 
"conventional" terms as given in e.g. [45]. Neither of 
these methods seems easy to be applied. But, without 
having fixed the value of /o at a sufficient precision, it 
is impossible to improve the accuracy of the prediction 
of the pion-deuteron scattering length. 

vi) In our opinion, it is feasible to estimate /g by using the 
resonance saturation hypothesis (see appendix G). At 
present time, however, the parameters of the resonance 
Lagrangian are not known at a sufficient precision, and 
more effort is needed to pin them down accurately from 
the experimental data. 
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A Calculation of the energy shift of pionic 
deuterium by using non-relativistic 
Lagrangians 

In this appendix, we shall present the calculation of the 
energy of the levels of the pionic deuterium, characterized 
by the quantum numbers nlj. The calculations will be 
performed up to and including ncxt-to-leading order in 
isospin-breaking parameters a and nid — m„. Since similar 
calculations for other hadronic bound states have been 
already considered in the literature in great detail (see e.g. 
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Fig. 9. The contributions to the electromagnetic energy of 
the pionic deuterium; a) relativistic insertions; b) one Coulomb 
photon exchange; c) one transverse photon exchange (shaded 
squares denote pion and deuteron electromagnetic formfac- 
tors); d) vacuum polarization. 



j"(p,q)^ 



r(p,q)=- 



2AL 



(A3) 



where Q = p + q, and = 0.2859 fm^, i^ij^i = 
2^Mm = 0.857412^ are the deuteron quadrupolc and 
magnetic moments, respectively. Substituting these rela- 
tions in Eq. (A2), we finally obtain 



[23-29]), we shall not repeat them here and display only 
the most important results. For convenience, the energy 
of the generic level can be split into the "electromagnetic" 
and "strong" parts 



(Al) 



where stands for the non-relativistic Coulomb binding 
energy. 

The contributions to the "electromagnetic" part of the 
potential are depicted in Fig. 9. These include: the rela- 
tivistic insertions p'^/SM^, p^/8M| in the pion and deute- 
ron lines, respectively; the one (Coulomb and transverse) 
photon exchange between the pion and the deuteron; the 
vacuum polarization contribution. At the end of the day, 
the explicit expression for the potential at this order can 
be written as 

Kr(P,q)-(2.)^^^(q-p)^..(^ + ^) 



+J°,(p,q)j^/(p,q) 

+^r.(p,q)5^(^..-^^^-^ 

-f<S,,y™^(p,q), 



/(P,q) 



(A2) 



where k = p — q and a,b = 1,2,3 denote the polarization 
projection of the deuteron. Further, stands for the 

vacuum polarization contribution and J^^, j'^ where /i = 
0, a denote the electromagnetic formfactor of the deuteron 
and the pion, respectively. At the order of accuracy we are 
working, one may use the following expression for these 
formfactors 



J°.(P,q) = e|^..(l - I (rl) k^) + ^ (k^h - Ik'S^b) 

+ ^{Qakb - kaQb){l - Mm)| + O(^) , 
J^biP, q) = ^ {Q'^^-b + flAliS^h - kaS^)} + O(^) , 



<5Kr(P,q) = (2^)''5^(q-p)'5a; 



/ 4 4 

I P : P 



k2 \ 6 
1 



Sab {{rl) + (rl)) 



q2 _ {Q^ ^ ) j 



k2 



2 



kakb 77 k Sab 



1 ^ , ^J^M 



(Qakb - kaQb) 



+ (5afcy™^(p,q). 



(A4) 



The electromagnetic energy shift at ncxt-to-lcading order 
is given by 



d^p d^q 

X - p)lp) Yil^_^^ {v)xl{o)Ki{\p\) 

X <5Kr(p,q)5Z^„;(|q|)x6(p)>l(._p)(q), (A5) 

where Y/^ and (jm|Z(m — <t)1<t) arc the spherical functions 
and the Clcbsch-Gordan coefficients, respectively and 

(AS) 

Note that the energy shift (A5) docs not depend on the 
magnetic quantum number v. Substituting here the ex- 
pression (A4), and calculating the integral, we finally ob- 
tain 



.pem_ Mj + MU apA f 4n 
"'^ MljMl \ n ) \l + l/2 
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6n 



10 



1 + 1/2 



AE%+AE% 



(A7) 



where for I ^ 



ApQ ( 30(2^-2)! 
- V (2^ + 3)! 



^ 1 J 
1 I 2 



(A8) 



4, ,3 



zli;„^f, = (-l)('+^)^2x/6 



1 - 
2MJ 



I 1 Z 
1 J 1 



(A9) 



v/Z(/ + l)(2/ + l) ' 

and for I = 0, AE^^^ = AE^J- = 0. The quantities in the 
braces denote the Wigner 6j symbols. Finally, the vacuum 
polarization contribution is given in Eq. (3) of Ref. [31]. 
We do not display it here. 
Taking / = 1, we have 



/ 



AE^, . 



3n3 



0.015 eV 



_ 1 

2 

J_ 

10 

/ 1 

.1 
2 



AE. 



3n3 



10 



Ma/ 



.7 = 0\ 

J = l 
.1 = 2) 

/ 1 



(AlO) 



2Mj MdM^ 



J = l 
.? = 2/ 



(All) 



As we see, the splitting of the energy levels is tiny. Note 
that in the averaged level energy, defined by Eq. (2.1), 

these contributions vanish: AE^^^ = AE:^Jj = 0. The re- 
maining contributions can be rearranged in order to allow 
the comparison with the results existing in the literature. 
In this manner, we finally arrive at the result given in Eqs. 
(2.2) and (2.3). 





( ' 


j 


= '\ 


0.027 eV 


1 




= 1 




2 


j 






j 


-2/ 



B Table of integrals 

After integrating over the 0th components A", /c", the quan- 
tity Ra in Eq. (3.22) can be rewritten as 



Ra = 



{2ttY {2Tr)d (72 + 12)(72_^k2) 

1 

72 + i[12+k2 + ^(l-k)2] • 



(Bl) 



Introducing Feynman parameters and carrying out the in- 
tegration over the momenta in the standard manner, one 
gets 

Ra^2m^{dl-2dl)- '^'^ ' ^ ^ 



(47r)« 



X / dxidx2dxs5{l — xi — X2 ~ x:i)[Ga] 



r, X3 X3 m 

G„ = xiXo H 1 

2 I 2 



(B2) 



where /i denotes the scale of dimensional regularization. 
Rcnormalizing Ra according to the MS prescription, for 
the finite part we get 



R. 



327r3 



Xa lnjj^ + Ja]+Oid- 3) 



Xa=- J dxidX2dX36{l - Xi- X2- X3)[GaY 

8M„ . m 



771 



arcsm ■ 



-1.2569- •• , 



m + Mj; 

Ja^^ J dxidX2dX3S{l - Xi-X2~ X3)[Ga]^i hj[Ga] 

(B3) 



= -1.1334- • • . 
In the same manner, one obtains 



(47r) 



X / dx{l — x) 

Jo \2hn 



1-^1 ,(B., 



where = mM-^im + M-,:)^^ is the reduced mass in 
the TriV system. After the renormalization one may evalu- 
ate the remaining integrals analytically. As the result, one 
obtains 



R 



^,,^ mHdi + 2di) 
327r3 



Xfc In ^ + J,, ) + 0{d - 3) , 



Xb-' 



:32k 5^1^= 0.5098 



Jb = 16k 2 1 ^^(1 — 2k) arcsin>/K 
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■ Vl - K (2 + ln(l - k)- In 4k) 



-0.3731 



2m 



(B5) 



The contribution Rc can be evaluated in a similar manner. 
We do not display this calculation here, since this contri- 
bution is proportional to and drops anyway. The same 
method can be applied to the calculation of the diagrams 
contributing to the threshold scattering amplitude A{fi). 



C Resonance saturation for the LEC /o 

In order to get an estimate for the LEC /g, we have cal- 
culated the tree amplitude of the process nNN ttNN 
in the theory with the explicit A, A''*(1440), scalar and 
vector mesons in the limit when masses of all above res- 
onances become very large. At threshold, the A and the 
vector meson do not contribute^. In Fig. 10 we display an 
example of a diagram that does not vanish at threshold. 
One has 4 such diagrams (insertion of iV*(1440) in each 
nucleon leg), as well as the diagrams which are obtained 
by a permutation of the outgoing nucleon legs. 

The Lagrangian describing the interactions between 
nucleons, pions, scalar mesons and the iV*(1440) is given 
by (see also [47]) 



C* = gsNNSNN + [g,NN'SN*N + h.c 



+ [clN*x+N~ 



idf,d^N*)u''u''N + h.c.], (CI) 



where x+ = M^{2 - tv'^/F'^ + ■ ■ ■) and = ~Td^n/F + 
• • •, and m* denotes the mass of the A''* (1440). The con- 
tribution to the scattering amplitude at threshold in the 
pertinent spin-isospin channel from the diagram in Fig. 10 
is given by 



-(2m) 



2 2M^(c* - C2)gsNNgsNN' 



X (xIxiX2X2 +xIx2X2Xi) , 



(C2) 



where Ms stands for the mass of the scalar meson. From 
the above expression, one may read off the contribution 
to the LEC 

■'° 2F^mAatiF ) ' 



■m*Ml 



(C3) 



where, according to the natural size arguments (see 
Eq. (4.14)), the quantity ^ must be of order 1. 

In the numerical estimates, we use the scalar meson 
mass Ms = 550 MeV (as it was also done in Ref. [47]). Fur- 
ther, gsNN and gsNN* denote the values of the pertinent 

^ Vector mesons are described by antisymmetric tensor fields. 



Fig. 10. Tree-level contribution to the scattering amplitude 
ttNN ttNN with the exchange of the Af*(1440) and the 
scalar meson. 



formfactors at zero momentum transfer. With the choice 
of monopolc formfactor Fsiq"^) = (vl| — M|)/(yl| — q"^) 
with As = 1700 MeV for both SNN and SNN* vertices, 
we deduce gsNN = 7.57 and gsNN* = 3.66 from the mass- 
shell values given e.g. in Ref. [48]. Further, in the same pa- 
per, one finds two different values for the constants c*, C2 



c* = -7.27 GeV"^ , c* = 



[Set 1] , 



-12.7 GeV" 



1.98 GcV"^ [Set 2] . 



(C4) 



Substituting these values, we get ^ = 21 [Set 1] and ^ = 
43 [Set 2] that obviously contradicts to our natural-size es- 
timate. Note that the values for c|,C2 given in Ref. [48] 
also do not agree with the result of the fit c-l + C2 = 
(— 1.56± 3.35) GeV~^ from Ref. [42]. However, in order 
to use the latter fit in our estimate, needs one more re- 
lation between cl and c^. We get this relation, assuming 
that these constants are saturated by the scalar meson 
exchange. The Lagrangian describing the interaction of 
pions with the scalar field with total isospin 7 = can be 
taken from Ref. [49] 



CdSTT(UnU^) 



^5Tr(x^ 



(C5) 



from the saturation hypothesis one gets c|/c2 = Cm/cd = 
4.2/3.2 (in the last equation, we have used numerical val- 
ues from Ref. [49]). Together with the value of cl + 
from Ref. [42] this yields the estimate ^ = 0.6 ± 1.3, which 
is in a reasonable agreement with the natural-size conjec- 
ture. Note that, albeit the resonance saturation method 
at present does not lead to an improved accuracy as com- 
pared to the dimensional analysis, the above discussion 
demonstrates the feasibility of estimating /q and eventu- 
ally /o through resonance saturation. What is needed to 
this end is to determine the LECs d{ , more accurately 
from the fit to A^*(1440) Ntttt and, in particular, to 
resolve the large discrepancies reported in Refs. [42,48]. 
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Erratum 

In the above paper we have studied the pion-deuteron 
scattering length within low-energy effective theories of 
QCD. In particular, in the section 4.3 of this article we give 
an estimate of the uncertainty in this scattering length, 
coming from the presence of the 6-particle low-energy con- 
stant (LEG) - referred to as /q. Unfortunately, in this 
estimate, which has yielded AaT^d/oL-Kd — 12%, an incor- 
rect value for the deuteron wave function at the origin 
in the Weinberg approach '^o(^) = [0.487; 0.434] fm~^/^ 
with the cutoff mass from the interval A — [450; 650] has 
been used (normalization error) . Using the corrected value 
#o(yl) = [0.137; 0.122] fm"^/^ downsizes the error (only 
due to the above-mentioned source) to a much more com- 
fortable AaTrd/a-wd — 1%, which, in addition, fits better 
to the a priori estimates, based on the Weinberg power 
counting. Thus, the results of the high-precision calcula- 
tions of Ref. [20] have been independently confirmed. 

This result has far-reaching consequences. In partic- 
ular, we have to re-think the equivalence between the 
Weinberg [20] and HPEFT [18, 19] approaches to the pion- 
deuteron scattering at low energy. The arguments remain 
the same, but the big change in the estimated uncertainty 
leads to the conclusion that now this equivalence is real- 
ized in a slightly different manner. As before, the success of 
the modified power counting [20] unambiguously indicates 
that HPEFT is a physically equivalent tool for describing 
pion-deuteron interactions near threshold. Furthermore, 
the sole input which HPEFT imports from the Weinberg 
approach are the values of LECs, determined through the 
matching procedure. As it is discussed in the above pa- 
per, these LECs in general acquire contributions from two 
different momentum regions: Af^ < p < A and p > A. 
Whereas the former can in principle be expressed through 
the wave functions, etc - and is thus calculable in terms 
of the known parameters of the theory, the latter repre- 
sents a genuinely high-energy contribution, which at the 
present stage can only be included in the error estimate. 
What follows from our corrected calculations is that this 
high-energy contribution is small and docs not preclude 
a high-precision determination of the pion-nucleon scat- 
tering lengths from the combined analysis of the pionic 
hydrogen and pionic deuterium data by Pionic Hydrogen 
collaboration at PSI [1,2]. We wish to emphasize once 
more that the strong predictive power is a result of a sub- 
tle balance between HPEFT, which at the end allows for 
constructing multiple scattering series in terms of observ- 
ables (scattering lengths, etc) and the Weinberg approach, 
which enables one to evaluate the LECs of HPEFT with 
a high precision. 

Several remarks are in order. First, we wish to mention 
that our error estimate, which is carried out by using di- 
mensional arguments or resonance saturation hypothesis, 
is complementary to the study of the scale dependence of 
the calculated pion-deuteron scattering length. Obviously, 
the investigations should proceed from both ends. As men- 
tioned in the above paper, a weak scale dependence would 
indicate that the uncertainty might be small and not that 



it must be small. Finally, the fact that the short-range in- 
teractions do not introduce large uncertainty in the calcu- 
lated value of the pion-deuteron scattering length, might 
stimulate further activity in calculating various small con- 
tributions, with an aim to determine this scattering length 
as precise as possible. Using the above-mentioned equiva- 
lence between the Weinberg and HPEFT approaches could 
enable one to carry out these calculations in a simpler set- 
ting. 



